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In this paper we study backward stochastic difTerential equations 
with general terminal value and general random generator. In par- 
ticular, we do not require the terminal value be given by a forward 
diffusion equation. The randomness of the generator does not need 
to be from a forward equation, either. Motivated from applications 
to numerical simulations, first we obtain the L*'-H61der continuity 
of the solution. Then we construct several numerical approximation 
schemes for backward stochastic differential equations and obtain the 
rate of convergence of the schemes based on the obtained L*'-H61der 
continuity results. The main tool is the Malliavin calculus. 



1. Introduction. The backward stochastic differential equation (BSDE, 
for short) we shall consider in this paper takes the following form: 



(1.1) Yt = C+ / f{r,Yr,Zr)dr - / ZrdWr, 0<t<T, 



where W = {Wt}o<t<T is a standard Brownian motion, ^ is the given ter- 
minal value and / is the given (random) generator. To solve this equation 
is to find a pair of adapted processes Y = {Yt}Q<t<T and Z = {Zt}o<t<T 
satisfying the above equation (1.1). 

Linear backward stochastic differential equations were first studied by 
Bismut [3] in an attempt to solve some optimal stochastic control problem 
through the method of maximum principle. The general nonlinear backward 
stochastic differential equations were first studied by Pardoux and Peng [15]. 
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Since then there have been extensive studies of this equation. We refer to 
the review paper by El Karoui, Peng and Quenez [7], to the books of El 
Karoui and Mazliak [6] and of Ma and Yong [12] and the references therein 
for more comprehensive presentation of the theory. 

A current important topic in the applications of BSDEs is the numerical 
approximation schemes. In most work on numerical simulations, a certain 
forward stochastic differential equation of the following form: 



is needed. Usually it is assumed that the generator / in (1.1) depends on Xy. 
at the time r: f{r,Yr,Zr) = f{r, Xr,Yr, Zr), where /(r, x, y, z) is a determin- 
istic function of (r, x, y, z), and / is global Lipschitz in (x, y, z). If in addition 
the terminal value ^ is of the form ^ = H^Xt), where ^ is a deterministic 
function, a so-called four-step numerical scheme has been developed by Ma, 
Protter and Yong in [11]. A basic ingredient in this paper is that the so- 
lution {Yt}o<t<T to the BSDE is of the form Yt = u{t,Xt), where u{t,x) is 
determined by a quasi-linear partial differential equation of parabolic type. 
Recently, Bouchard and Touzi [4] propose a Monte-Carlo approach which 
may be more suitable for high-dimensional problems. Again in this forward- 
backward setting, if the generator / has a quadratic growth in Z, a numerical 
approximation is developed by Imkeller and Dos Reis [9] in which a trunca- 
tion procedure is applied. 

In the case where the terminal value ^ is a functional of the path of the 
forward diffusion X, namely, = g{X.), different approaches to construct 
numerical methods have been proposed. We refer to Bally [1] for a scheme 
with a random time partition. In the work by Zhang [16], the L^-regularity 
of Z is obtained, which allows one to use deterministic time partitions as well 
as to obtain the rate estimate (see Bender and Denk [2], Gobet, Lemor and 
Warin [8] and Zhang [16] for different algorithms). We should also mention 
the works by Briand, Delyon and Memin [5] and Ma et al. [10], where the 
Brownian motion is replaced by a scaled random walk. 

The purpose of the present paper is to construct numerical schemes for 
the general BSDE (1.1), without assuming any particular form for the termi- 
nal value ^ and generator /. This means that ^ can be an arbitrary random 
variable, and f{r,y,z) can be an arbitrary J>-measurable random variable 
(see Assumption 2.2 in Section 2 for precise conditions on and /). The 
natural tool that we shall use is the Malliavin calculus. We emphasize that 
the main difficulty in constructing a numerical scheme for BSDEs is usually 
the approximation of the process Z. It is necessary to obtain some regular- 
ity properties for the trajectories of this process Z. The Malliavin calculus 
turns out to be a suitable tool to handle these problems because the random 



(1.2) 
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variable Zi can be expressed in terms of the trace of the Mahiavin derivative 
of Yf, namely, Zt = DfYf. This relationship was proved in the paper by El 
Karoui, Peng and Quenez [7] and was used by these authors to obtain esti- 
mates for the moments of Zf. We shall further exploit this identity to obtain 
the L^-Holder continuity of the process Z, which is the critical ingredient 
for the rate estimate of our numerical schemes. 

Our first numerical scheme was inspired by the paper of Zhang [16], 
where the author considers a class of BSDEs whose terminal value takes 
the form g{X.), where X is a forward diffusion of the form (1.2), and g 
satisfies a Lipschitz condition with respect to the L°° or norms (simi- 
lar assumptions for /). The discretization scheme is based on the regular- 
ity of the process Z in the mean square sense; that is, for any partition 
TT = {0 = to < ti < ■ ■ ■ < tn = T} , one obtains 

(1.3) [\Zt-Ztf + \Zt-Zt^^,\^]dt<K\Tr\, 

where |7r| = maxo<j<n-i(ti+i — ti), and if is a constant independent of the 
partition vr. 

We consider the case of a general terminal value ^ which is twice differen- 
tiable in the sense of Malliavin calculus, and the first and second derivatives 
satisfy some integr ability conditions; we also made similar assumptions for 
the generator / (see Assumption 2.2 in Section 2 for details). In this sense 
our framework extends that of [13] and is also natural. In this framework, 
we are able to obtain an estimate of the form 

(1.4) E\Zt- Zs\P <K\t-s\P/^, 

where if is a constant independent of s and t. Clearly, (1.4) with p = 2 
implies (1.3). Moreover, (1.4) implies the existence of a 7-IIolder continuous 
version of the process Z for any 7 < ^ ~ |- Notice that, up to now the 
path regularity of Z has been studied only when the terminal value and the 
generator are functional of a forward diffusion. 

After establishing the regularity of Z, we consider different types of nu- 
merical schemes. First we analyze a scheme similar to the one proposed 
in [16] [see (3.2)]. In this case we obtain a rate of convergence of the follow- 
ing type: 

E sup \Yt-Yt''\'^+ [ E\Zt- Z]^\'^dt<K{\^T\+E\^-C\'^)■ 
o<t<T Jo 

Notice that this result is stronger than that in [16] which can be stated as 
(when ^'^ = ^) 

sup E\Yt-Yt''\'^ + [ E\Zt- Z]^\'^dt<K\7r\. 

0<t<T Jo 
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We also propose and study an "implicit" numerical scheme [see (4.1) in 
Section 4 for the details] . For this scheme we obtain a much better result on 
the rate of convergence, 

/ i-T \p/2 
E sup \Yt-Yt''\P + Ei / iZt-Zfpdt <K{\Tr\P/^ +E\^-Cn, 

0<t<T \Jo J 

where p> 1 depends on the assumptions imposed on the terminal value and 
the coefficients. 

In both schemes, the integral of the process Z is used in each iteration, 
and for this reason they are not completely discrete schemes. In order to 
implement the scheme on computers, one must replace an integral of the 
form Jf.'^^ Zg ds by discrete sums, and then the convergence of the obtained 
scheme is hardly guaranteed. To avoid this discretization we propose a truly 
discrete numerical scheme using our representation of Zt as the trace of the 
Malliavin derivative of Yt (see Section 5 for details). For this new scheme, 
we obtain a rate of convergence result of the form 

E max {\Yu - Y,^ + \Zu - Z^ < K\7r\P/^-^ 

for any e > 0. In fact, we have a slightly better rate of convergence (see 
Theorem 5.2), 

E max {\Yt^ - y^^l^ + \Zt^ - Zl\^} < K|^|P/2-p/(2iog(i/k|)) ^og J- 

0<i<n » 1 » I \^ 

However, this type of result on the rate of convergence applies only to some 
classes of BSDEs, and thus this scheme remains to be further investigated. 

In the computer realization of our schemes or any other schemes, an ex- 
tremely important procedure is to compute the conditional expectation of 
form E(y|J^t. ). In this paper we shall not discuss this issue but only mention 
the papers [2, 4] and [8]. 

The paper is organized as follows. In Section 2 we obtain a representa- 
tion of the martingale integrand Z in terms of the trace of the Malliavin 
derivative of Y, and then we get the L^-Holder continuity of Z by using 
this representation. The conditions that we assume on the terminal value ^ 
and the generator / are also specified in this section. Some examples of ap- 
plication are presented to explain the validity of the conditions. Section 3 
is devoted to the analysis of the approximation scheme similar to the one 
introduced in [16]. Under some differentiability and integrability conditions 
in the sense of Malliavin calculus on ^ and the nonlinear coefficient /, we 
establish a better rate of convergence for this scheme. In Section 4, we in- 
troduce an "implicit" scheme and obtain the rate of convergence in the U' 
norm. A completely discrete scheme is proposed and analyzed in Section 5. 

Throughout the paper for simplicity we consider only scalar BSDEs. The 
results obtained in this paper can be easily extended to multi-dimensional 
BSDEs. 
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2. The Malliavin calculus for BSDEs. 

2.1. Notations and preliminaries. Let W = {Wt}o<t<T be a one-dimen- 
sional standard Brownian motion defined on some complete filtered proba- 
bility space {Q,J^,P, {Tt}o<t<T)- We assume that {Tt}o<t<T is the filtration 
generated by the Brownian motion and the P-null sets, and T = Tt- We 
denote by V the progressive u-field on the product space [0, T] x 0. 

For any p > 1 we consider the following classes of processes: 

• M^'P, for any p>2, denotes the class of square integrable random variables 
F with a stochastic integral representation of the form 

F = EF+ [ utdWt, 
Jo 

where n is a progressively measurable process satisfying supo<t<TE|nt|^ < 
oo. 

• H^{[0,T]) denotes the Banach space of all progressively measurable pro- 
cesses f: {[0,T] xn,!^)^ {R,B) with norm 

/ / /-T xp/2n 1/p 

\\ip\\HP = [^ij Inl'^dtj ] <oo. 

• Sjr{[0,T]) denotes the Banach space of all the ROLL (right continuous 
with left limits) adapted processes ip:{[0,T] x il.,V) (K, S) with norm 

/ N i/p 

||(^||sp = (E sup liptl''') < oo. 

Next, we present some preliminaries on Malliavin calculus, and we refer 
the reader to the book by Nualart [14] for more details. 

Let H = L2([0,T]) be the separable Hilbert space of all square integrable 
real- valued functions on the interval [0,T] with scalar product denoted by 
(■)')h- The norm of an element /i S H will be denoted by ||/i||h- For any 
/i e H we put W{h) = /(f h{t) dWt. 

We denote by C^(M") the set of all infinitely continuously differentiable 
functions g : — ?> M such that g and all of its partial derivatives have poly- 
nomial growth. We make use of the notation dig = ^ whenever g € C"^(M"'). 

Let S denote the class of smooth random variables such that a random 
variable F € 5 has the form 

(2.1) F = g{W{hi),...,W{K)), 

where g belongs to C^(M'"), /ii, . . . , /i„ are in H and n > 1. 

The Malliavin derivative of a smooth random variable F of the form (2.1) 
is the H-valued random variable given by 

n 

DF = Y, ^^9iW{hl), W{hn))hi. 
i=l 
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For any p > 1 we will denote the domain of D in L^{il) by 0"*^'^, meaning 
that B^'^ is the closure of the class of smooth random variables S with 
respect to the norm 

||Fi|i,p = (E|F|P + E||Z)F||fj)i/P. 

We can define the iteration of the operator D in such a way that for a smooth 
random variable F, the iterated derivative D^F is a random variable with 
values in H®'^. Then for every p > 1 and any natural number > 1 we 
introduce the seminorm on S defined by 

i/p 

\F\ ' ^ . ^ \ 



\k,p ■ 



/ k \ 

i^E\F\^+j2nD^Fr^^^j 



We will denote by B'^'*' the completion of the family of smooth random 
variables S with respect to the norm |[ • \\k,p- 

Let /i be the Lebesgue measure on [0,T]. For any k>l and F G B'^'*', the 
derivative 

D''F = {Dt^_,^F,Ue[0,T],i = l,...,k} 

is a measurable function on the product space [0,T]'^ x 0, which is defined 
a.e. with respect to the measure fi^ x P. 

We use Lq^ to denote the set of real-valued progressively measurable 
processes u = {ut}o<t<T such that: 

(i) For almost ah t € [0,T],ut € B^'P. 

(ii) E((/o^ \ut\^dt)P/^ + (/o^/o^ \DeUt\'dedt)P/^) < oo. 

Notice that we can choose a progressively measurable version of the H- valued 
process {Dut}o<t<T- 

2.2. Estimates on the solutions of BSDEs. The generator / in the BSDE 
(1.1) is a measurable function /: ([0,r] x Q x R x R,V ^ B ^ B) ^ (R,B), 
and the terminal value is an J-T^-measurable random variable. 



Definition 2.1. A solution to the BSDE (1.1) is a pair of progressively 

[q \Zt\'^dt <oo, 



measurable processes {Y, Z) such that \Zt\'^ dt<oo, \f{t,Yt, Zt) \ dt <oo, 



a.s. and 

rT 



Yt = ^+ [ f{r,Yr,Zr)dr- [ ZrdWr, 0<t<T. 
Jt Jt 

The next lemma provides a useful estimate on the solution to the BSDE (1.1). 



MALLIAVIN CALCULUS, NUMERICAL SOLUTION OF BSDE 



7 



Lemma 2.2. Fix q>2. Suppose that ^ G Li{n), f{t, 0,0) e F|-([0,r]) 
and f is uniformly Lipschitz in {y,z); namely, there exists a positive num- 
ber L such that fi x P a.e. 

\f{t,yi,zi) - f{t,y2,Z2)\ < L{\yi - y2\ + \zi - Z2I) 

for all yi,y2 G M and zi,Z2 G K. Then there exists a unique solution pair 
{Y,Z) G 5^([0,T]) X HmO,T]) to (1.1). Moreover, we have the following 
estimate for the solution: 

/ rT \ q/2 

E sup l^tl^ + E / \Zt\'^dt 



0<t<T \Jo 

(2.2) 



<k(k\C\'^ + e(^J^ \f{t,0,0)\^dt^ 



where K is a constant depending only on L, q and T. 

Proof. The proof of the existence and uniqueness of the solution (Y, Z) 
can be found in [7], Theorem 5.1, with the local martingale M = 0, since the 
filtration here is the filtration generated by the Brownian motion W. Esti- 
mate (2.2) can be easily obtained from Proposition 5.1 in [7] with (/^,^^) = 
(/,Oand(/2,a = (0,0). □ 

As we will see later, for a given BSDE the process Z will be expressed in 
terms of the Malliavin derivative of the solution Y, which will satisfy a linear 
BSDE with random coefficients. To study the properties of Z we need to 
analyze a class of linear BSDEs. 

Let {at}o<t<T and {/3t}o<t<T be two progressively measurable processes. 
We will make use of the following integr ability conditions: 

Assumption 2.1. 
(HI) For any A > 0, 



Cx :=Eexp^A^ 



T 

at\ + I3t)dt ] < 00. 



(H2) For any p > 1, 

Kp:= sup E{\at\P + \(3tn < 00. 

0<t<T 

Under condition (HI), we denote by {pf}o<t<T the solution of the linear 
stochastic differential equation 

' dpt = atptdt + /3tptdWt, 0<t<T, 

/Oo = l. 

The following theorem is a critical tool for the proof of the main theorem 
in this section, and it has also its own interest. 



(2.3) 
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Theorem 2.3. Let q> p>2 and let ^ e L^iQ) and f £ H'^{[0,T]). As- 
sume that {aj}o<t<T CLnd {(3t}o<t<T o-'re two progressively measurable pro- 
cesses satisfying conditions (HI) and (H2) in Assumption 2.1. Suppose that 

the random variables (,pT and Jq ptftdt belong to M^'"^, where {pt}o<t<T is 
the solution to (2.3). Then the following linear BSDE, 

(2.4) Yt = i+ [ [arYr + l3rZr + fr]dr - I ZrdWr, 0<t<T, 

Jt Jt 

has a unique solution pair (Y,Z), and there is a constant K > such that 

(2.5) E\Yt-Ys\P <K\t-s\P/^ 
for all s,te [0,T]. 

We need the following lemma to prove the above result. 

Lemma 2.4. Let {at}o<t<T and {(3t}o<t<T be two progressively measur- 
able processes satisfying condition (HI) in Assumption 2.1, and {pt}Q<t<T 
be the solution of (2.3). Then, for any r € M i(;e have 

(2.6) E sup yo[<oo. 

0<i<T 

Proof. Let t G [0, T]. The solution to (2.3) can be written as 

pt = expj^ {as - y) + _^ 
For any real number r, we have 

E sup /o[=E sup exp| / r(as-^^^ds + r I fisdWg 
o<t<T o<t<T [Jo \ 2 / Jo 

E^exp||r|^ |Qs|ds+ ^(|r| +r2)^ 



< 



0<t<T 



t ^2 ft 



2 



X sup exp<! r / PsdWs — — I P^ds 



Then, fixing any p> 1 and using Holder's inequality, we obtain 
(2.7) E sup pI<c(w. sup expjrp / PsdWs-^ f P^ds 

0<t<T \ 0<t<T [Jo 2 Jo 

where 

C= ^Eexp|g|r|^ ds + |(|r| + r^) ^ /jfdsj^ 



i/p 



and ^ + i = l. 

p q 
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Set Mt = exp{r /3s dWs — ^ Jq (3l ds}. Then {Mt}o<t<T is a martingale 
due to (HI). We can rewrite (2.7) into 

i/p 



(2.8) 



E sup p[<C E sup Mf 

0<t<T ^ 0<t<T 



By Doob's maximal inequality, we have 



(2.9) 



E sup Mf<CpEM^ 

0<t<T 



for some constant Cp > depending only on p. Finally, choosing any 7 > 1, 

7 + A 



A > 1 such that - + 4 = 1 and applying again the Holder inequality yield 



EM^ = E (^exp I rp dWs - 



T 



X exp 



7P- 1 2 
— - — pr 



Plds 



Plds 



< ( Eexp<j rpj I PsdWs — - j'^p'^r'^ 



X Eexp 



KlP - 1) 2 



(31 ds 



l/A 



1/7 



Eexp 



KlP - 1) 2 
pr 



f3lds 



l/A 



< 00. 



Combining this inequality with (2.8) and (2.9) we complete the proof. □ 

Proof of Theorem 2.3. The existence and uniqueness is well known. 
We are going to prove (2.5). Let t € [0,T]. Denote 74 = p~[^ , where {pt}f)<t<T 
is the solution to (2.3). Then {7t}o<t<r satisfies the following linear stochas- 
tic differential equation: 

dit = {-at + /3|)7t dt - Pat dWt, 0<t<T, 
70 = 1- 

For any < s <t <T and any positive number r > 1, we have, using (H2), 
the Holder inequality, the Burkholder-Davis-Gundy inequality and Lem- 
ma 2.4 applied to the process {'yt}o<t<T, 



E|7t - 7s 



E 



(2.10) 



< 2 



E 



-Ou + Plhu du - 



-ttu+Plhu du 



Pulu dWu 

r 

+aE 



< c{t - sy/^, 
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where Cr is a constant depending only on r, and C is a constant depending 
on T, r and the constants appearing in conditions (HI) and (H2). 
From (2.3), (2.4) and by Ito's formula, we obtain 

d{Ytpt) = -ptft dt + {ptPtYt + ptZt) dWt. 

As a consequence. 



(2.11) Yt = Pt^¥.(^ipT + prfrdr\F^=¥.(^ipt,T + Pt,rfrdr 

where we write pt^r = Pt^ Pt = ItPr for any < t < r < T. 
Now, fix < s < t < T. We have 



Ely, - yjp = E 



T 



E iPt,T+ / Pt,rfrdr Tt] -E( ^Ps,T+ / Ps,rfrdr 



< 2^- 



E|E(Car|-Ft)-E(ep.,T|-F.)|^ 



+ E 



E 



Pt,rfrdr 



E 



Ps,rfr dv 



PI 



First we estimate I\. We have 
/i =E|E(^p,,T|^t) -IE(ep.,T|^.)|^ 

= E|E(C,9t,T|-Ft) - E(^p,,T|-Ft) + E(ep.,r|-Ft) - E(ep.,r|-F.)r 

< 2P-i[E|E(e/>t,T|-^t) -IE(ePs,T|-^t)r + IE|E(Cps,T|-^t) -IE(CPs,T|^.)r] 

< 2P-nE|e(aT - Ps,T)r + IE|E(C/5,,T|-^t) - mPs,T\:Fs)n 
= 2P-\h + h). 

Using the Holder inequality. Lemma 2.4 and the estimate (2.10) with r = 
the term Is can be estimated as follows: 

/3 < (Eicr)p/'^(E|/>,,T - 

< (E|^|'?)f/'?(E|7t - ^^|W(<?-p))(5-p)/{2g)^E^2pg/(g-p)^(g-p)/{2g) 

<C7|i-s|P/2, 

where C is a constant depending only on p,q,T, E|^|'? and the constants 
appearing in conditions (HI) and (H2). 

In order to estimate the term I4 we will make use of the condition S^pT € 
M^''^. This condition implies that 

^PT = Ei^PT)+ [ UrdWr, 

Jo 
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where u is a progressively measurable process satisfying supo<«rIE|tit|'' < oo. 
Therefore, by the Burkholder-Davis-Gundy inequality, we have 

EmpT\Tt)-E{MJ's)\' 



E 



Ur dWr 



u^. dr 



q/2 



<Cq{t-sY''^ sup E|nt|«. 

Q<t<T 

As a consequence, from the definition of I4 we have 

h = E\^MiPT\:Ft) - E{ipT\FsW 

<C\t-s\P/\ 

where C is a constant depending on p,q,T,supQ^f^rpE\ut\'^ < 00 and the 
constants appearing in conditions (HI) and (H2). 
The term I2 can be decomposed as follows: 



T 



E / pt,rfrdr J='t] -El / Ps,rfrdr 



T 



E 



Pt,rfrdr Tt^ -E(^j^ ps^rfrdr Tt^ 



Tt] - El / Ps,rfrdr 



Tt] - El / Ps,rfrdr 



= 3P-\h + h + l7)- 

Let us first estimate the term I^. Suppose that p <p' < q. Then, using (2.10) 
and the Holder inequality, we can write 



h = E 



< E 



Pt,rfrdr 

T 

{Pt,r - Ps,r)frdr 



Tt]-E 



Ps,rfr dr 



IE It* -7s 



< 



{E|7t - 7,|PP'/{p'-p)}(p'-p)/p'i]E 



T 



Prf r dr 



Prfrdr 
P' ^ P/p' 
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'V 



.T \p'q/{-2{q-p'))^p{q-p')/{p'q) 

<C7|t_s|p/2j]E / pldr 



rT \q/2sp/q 

X <jE( / /2drj I 



<C\t-sf^\\f\\ 



mi 



where C is a constant depending on q, T and the constants appearing 
in conditions (HI) and (H2). 

Now we estimate Jg. Suppose that p <p' < q. We have, as in the estimate 
of the term I5 , 



<E 



T 

E( / Ps,rfrdr 
P 

Ps,rfr dr 



T 

-Ft ) -E( / Ps,rfrdr 

p 

Prfrdr 



^[Pl" 

< {E/>7Pp'/(p'-p)}(p'-p)/p'|e J' Prfrdr 

Prfrdr 
|P/2| 



p' ~\ p/p' 



p' ^ p/p' 



<C|*-.|-1E sup f"'""'"||/ll'„. 

0<t<T ^ 

d\t-s\p/\ 



where C is a constant depending on p,p' , q, T and the constants appearing 
in conditions (HI) and (H2). 



The fact that Prfrdr belongs to M^'"^ imphes that 



Prfrdr = 'E / Prfrdr + / VrdWr, 

Jo Jo 

where {vt}o<t<T is a progressively measurable process satisfying 

sup E|fj I'' < 00. 
o<t<r 

Then, by the Burkholder-Davis-Gundy inequality we have 



E 



E / Prfrdr 



-E / Prfrdr 



T 



E / Prfrdr 



T 



Tt]-m / Prfrdr 



Vr dWr 



<Cg{t-sy/^ sup Elvtl". 

0<t<T 
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Finally, we estimate as follows: 

rT 



I- 



7 



E 



E 



E 



Ps,rfr dv 
T 



E 



T 



E 



Prfr dr 



E 



Ps,rfr dv 
T 



Prfrdr 



(2.12) 



rT 



X <^ E 



E 



Prfr dr 



■E 



Prfrdr 



<C<^E 



E 



T 



Prfrdr 



■E 



Prfrdr 



<i-\ pN 



<C|t-s|P/2, 



where C is a constant depending on p, T, supo<t<TE|'yt|'' and the con- 
stants appearing in conditions (HI) and (H2). 
As a consequence, we obtain for all s,t € [0, T] 

E\Yt-Ys\P <K\t-s\P/^, 

where K is a constant independent of s and t. □ 

2.3. The Malliavin calculus for BSDEs. We return to the study of (1.1). 
The main assumptions we make on the terminal value ^ and generator / are 
the following: 

Assumption 2.2. Fix2<p<|. 
(i) i £ D^.g^ i-j^g^g gxists L > 0, such that for ah 6,6' G [0,r], 
(2.13) E\DeC- De'C\P <L\e-6'\P/'^, 



(2.14) 

and 
(2.15) 



sup El^e^l" < oo 

0<6KT 



sup sup E,\DuDg^\'^ < OO. 
0<6KT0<«<T 



(ii) The generator f{t, y, z) has continuous and uniformly bounded first- 
and second-order partial derivatives with respect to y and z, and /(■,0,0) € 
HU[0,T]). 

(iii) Assume that ^ and / satisfy the above conditions (i) and (ii). Let 
(Y, Z) be the unique solution to (1.1) with terminal value and generator /. 
For each {y,z) G M x M, f{-,y,z), dyf{-,y,z) and dzf{-,y,z) belong to La'^, 
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and the Malliavin derivatives Df{-,y,z), Ddyf{-,y,z) and Ddzf{-,y,z) sat- 
isfy 

/ rT \q/2 

(2.16) sup E / \Def{t,Yt,Zt)\^dt) < oo, 
o<e<r \Je / 

/ rT \<?/2 

(2.17) sup E / \Dedyf{t,Yt,Zt)\^dt) < oo, 
o<0<T \Je / 

/ rT \q/2 

(2.18) sup E / \Ded,f{t,YuZt)\''dt] < oo, 
0<6»<T \Je J 

and there exists L > such that for any t G (0, T] , and for any 0<9,9' <t<T 

f rT \p/2 

(2.19) EU \Dej(T,Yr,Zr)-De,f{T,Yr,Zr)fdT\ <L\e-e'\^l^. 

For each ^ G [0,r], and each pair of {y,z), Dof{-,y,z) G La"^ and it has 
continuous partial derivatives with respect to y,z, which are denoted by 
dyDgflt, y, z) and dzD^fit, y, z), and the Mahiavin derivative DuDgfit^ y, z) 
satisfies 

/ rT \g/2 

(2.20) sup sup E / \DuDef{t,Yt,Zt)\^dt] < oo. 
o<e<To<«<T \jeyu J 

The following property is easy to check and we omit the proof. 

Remark 2.5. Conditions (2.17) and (2.18) imply 

/ rT \g/2 
sup E / \dyDef{t,YuZt)\^dt\ < oo 
0<6»<T \Je J 

and 

/ rT \q/2 

sup E / \dzDgf{t,Yt,Zt)\^dt) <oo, 
o<e<T \Je / 

respectively. 

The following is the main result of this section. 
Theorem 2.6. Let Assumption 2.2 he satisfied. 

(a) There exists a unique solution pair {{Yt, Zt)}o<t<T to the BSDE (1.1), 
and Y,Z are in La*^. A version of the Malliavin derivatives {{DgYt, 
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DQZt)}Q<e,t<T of the solution pair satisfies the following linear BSDE: 



DeYt = Dei + j\dyf{r,Yr,Zr)DeYr 



(2.21) + d-J{r,Y,, Zr)DeZr + Def{r,Yr,Zr)] dr 
rT 

DgZrdWr, 0<0<t<T; 

(2.22) DeYt = 0, DeZt = 0, 0<t<e<T. 

Moreover, {DtYt}o<t<T defined by (2.21) gives a version of {Zt}o<t<T, namely, 
H X P a.e. 

(2.23) Zt = DtYt. 

(b) There exists a constant K > 0, such that, for all s,i G [0,r], 

(2.24) E\Zt- Z,\P <K\t- sf/"^. 

Proof. Part (a): The proof of the existence and uniqueness of the solu- 
tion (y, Z) , and Y,Z£ L^' is similar to that of Proposition 5.3 in [7] , and also 
the fact that {DeYt,DgZt) is given by (2.21) and (2.22). In Proposition 5.3 
in [7] the exponent q is equal to 4, and one assumes that \\D0f{-,Y, 

Z)\\'jj2 dO < oo, which is a consequence of (2.16) and the fact that Y,Z^ La^. 

Furthermore, from conditions (2.14) and (2.16) and the estimate in Lem- 
ma 2.2, we obtain 

(2.25) sup |e sup \DgYt\'^+E( [ {DeZtl"^ dt^ 1 < oo. 
o<e<T I e<t<T \Je J J 

Hence, by Proposition 1.5.5 in [14], Y and Z belong to La"^. 

Part (b): Let < s <t <T. In this proof, C > will be a constant inde- 
pendent of s and t, and may vary from line to line. 

By representation (2.23) we have 

(2.26) Zt-Z, = DtYt - DsYs = {D^Yt - D.Y^) + {D,Yt - DsY^). 

From Lemma 2.2 and equation (2.21) for 9 = s and 9' = t, respectively, we 
obtain, using conditions (2.13) and (2.19), 

/ j-T \p/2 

E\DtYt - DsYtf + eU \DtZr - D^Zr]^ drj 



<C 



(2.27) 



E\Dti-Dsi\ 



p 



/ rT \p/2 
+ E ( / \Dtf{r, Yr, Zr) - Dsf{r, Y^, Zr)f dr ] 

<C\t-s\P/\ 
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Denote a„ = dyf{u, Y^, Zy) and j3u = dzf{u, Yu, Z^) for all u G [0, T]. Then, 
by Assumption 2.2(ii), the processes a and /3 satisfy conditions (HI) and (H2) 
in Assumption 2.1, and from (2.21) we have for r G [s,T] 

D,Yr = Dsi+ f [auDsYu + PuDsZu + Dsf{u,Yu,Zu)]du- f D^ZudWu. 

J r J r 

Next, we are going to use Theorem 2.3 to estimate E,\DsYt — DsYs\p. Fix p' 
with p < p' < I (notice that p' < f is equivalent to < !)• From con- 
ditions (2.14) and (2.16), it is obvious that DsS, G ^^(0) C LP'(O) and 
Dsf{-,Y,Z) G H'^iiO^T]) C HP'{[0,T]) for any s G [0,T]. We are going to 
show that, for any s G [0, T], prDgi and /J" puDsf{u, Yu, Zu) du are elements 
in M'^'P' , where 

/Jr = exp|y (3udWu + j (^au - ^(3l^ du^ . 
For any < < r < T, let us compute 

Dgpr = Pr| y [dyzf{u,Yu,Zu)DgYu 

+ 5,,/(n, K,, Z„)Z)e^„ + DedJ{u, Y^, ZJ] dH^, 
+ dJ{e,Ye,Zg) 

rr 

+ / {9yyf{u, Yu, Zu) - dyj{u, Yu, Zu)Pu)DgYu du 
Je 

/•r 

+ / {dyzfiu, Yu, Zu) - dzzfiu, Yu, Zu)(3u)DgZu du 
Je 

+ j\Dedyf{u,Yu,Zu) - (3uDedzf{u,Yu,Zu))duy 

By the boundedness of the first- and second-order partial derivatives of / 
with respect to y and z, (2.17), (2.18), (2.25), Lemma 2.4, the Holder in- 
equality and the Burkholder-Davis-Gundy inequality, it is easy to show that 
for any p" < q, 

(2.28) sup E sup \Dg prf <oo. 

0<e<T e<r<T 

By the Clark-Ocone-Haussman formula, we have 

PTDs( = EipTD,C) + r nDeiPTDsOmdWe 
Jo 

= E{ptDsO+ [ E{DePTDsi + PTDeD,i\Fe)dWe 
Jo 

= E{ptD,0+ ruldWe 
Jo 



MALLIAVIN CALCULUS, NUMERICAL SOLUTION OF BSDE 



and 



/ prDsf{r,Yr,Zr)dr 

J s 

= E^ PrDJ{r,Yr,Zr)dr 

+ j\[De j%rDsf{r,Yr,Zr)dr\F^dWe 
= E^ PrDsf{r,Yr,Zr)dr 

+ j\(^£[DePrDsf{r, F,, Z,,) 

+ PrdyDJ{r,Yr,Zr)DeYr 
+ prdzDsf{r, Yr, Zr)DeZr 

+ PrDeDsfir,Yr,Zr)]dr 
= E [ prDJ{r,Yr,Zr)dr+ [ vIdWe. 

Js Jo 

We claim that supo<5i<r < oo and supQ<g<y E|'y||''' < oo. In fact 
E\u$\P' = E\E{DgpTDsC + prDeDSJ^e)]"' 

< 2P'~\E\DepTDs^f +E\pTDeDsCf) 
<2P'-\{E\DepTf'i/'^'i-P'^)^'^~^''^^%E\Ds^\Y'/i 

By (2.14), (2.15), (2.28) and Lemma 2.4, we have supo<5<rSupo<e<rE|ti 
oo. On the other hand, 



J'e dWe 



E\vf,\P = E 



E 



[ [DePrDsf{r,Yr,Zr 
J s 



+ PrdyDsf{r,Yr,Zr)DeYr 

+ prdzDsf{r,Yr, Zr)DeZr 
+ PrDeDsfir,Yr,Zr)]dr 
<4P'-i[Ji + J2 + J3 + J4], 
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Ji =E 



J2=IE 



J3 = E 



/ DePrDsf{r,Yr,Zr)dr 

J s 

[ prdyDsf{r,Yr,Zr)DeYrdr 

J s 

/ pr Dsf{r, Yr , Zr ) Dq Zr dr 

J s 



and 



For Ji, we have 



J4 = E 



PrDgDsf{r,Yr,Zr)dr 



Ji<e( sup \Dgprf [ Dsf{r,Yr,Zr)dr 

\e<r<T Js 

<(e sup |i;,p.PV(.-.')^^-^')/^ 



e<r<T 



X E 



r DJir,Yr,Zr)dr 

J s 



1\ p'/<l 



<r^''/2(E sup |z?,p.r'«/(«-p'))^''"''^^' 

^ e<r<T ' 

/ / /-T \g/2Np7g 
X (^E^ |D,/(r,y„Z,)|2drj j . 



For J2, we have 



J2<Ef sup \DeYr\^'{ sup / jSyL'./Cr.y^Z,)! drV 

Ve<r<T Vo<r<T Js / 



< E sup \DeYr\^ 

e<r<T 



v'li 



sup Pr [ \dyDsf{r,Yr,Zr)\dr 

V \0<r<T Js 



<(e sup \DeYrVY'"(^ sup /./(.-V)^ ^-^^'^/^ 

^ ^ 0<r<T 



e<r<T 



X (^e(^^ |9^Z),/(r,y,,Z,)|(iry^ 
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<T^'/'(e sup \DeY,.\A''^'(E sup ^-^^'^/^ 

^ e<r<T ' ^ 0<r<T 



E 



Using a similar techniques as before, we obtain that 
X (^E^ |a,Z),/(r,y,,Z,)pdrj j 



and 

(<?-p')/g 



J4<Tp'/2(e sup 

^ 0<r<T 



X (^E^ |L>eI?s/(r,i;,Z,,)NrJ J . 

By (2.16), (2.17)-(2.20), (2.28) and Lemma 2.4, we obtain that 

sup sup E|u||^ < oo. 

0<s<TQ<B<T 

Therefore, pxi and puDsf{u,Yu, Zu) du belong to M^'^'. 

Thus by Theorem 2.3 with p < p' , there is a constant C{s) > 0, such that 

E\DsYt-DsYs\P<C{s)\t-s\P/^ 

for all t G [SiT]. Furthermore, taking into account the proof of the esti- 
mates Ifc {k = 3,4, ... ,7) in the proof of Theorem 2.3, we can show that 
supo<s<T C'('S) =: C < oo. Thus we have 

(2.29) E\DsYt-DsYs\f<C\t-s\f^'^ 

for ah s,te [0,T]. Combining (2.29) with (2.26) and (2.27), we obtain that 
there is a constant K > independent of s and t, such that 

E\Zt-Zs\P<K\t-s\P/^ 

for aU s,tG [0,r]. □ 

Corollary 2.7. Under the assumptions in Theorem 2.2, let (Y,Z) € 
S^([0,T]) X i7^([0,T]) he the unique solution pair to (l-l)- If 
supo<t<T^I-^t|'^ < oo, then there exists a constant C, such that, for any 
s,te[Q,T], 

(2.30) E\Yt-Ys\'' <C\t- s\'^l'^. 
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Proof. Without loss of generality we assume 0<s<t<T. C>Ois 
a constant independent of s and t, which may vary from line to line. Since 

Ys = Yt+ f f{r,Yr,Zr)dr- f Zr-dWr, 

J s J s 

we have, by the Lipschitz condition on /, 



/ f{r,Yr,Zr)dr- j ZrdWr 

J s J s 



< 2''~M E 



ZrdWr 

\ 9/2 



f f{r,Yr,Zr)dr\ : 

J s 

<Cq[\t-s\''/^¥.(^j\f{r,Yr,Zr)\^dry' +E(^j\zr\^ drj^ ^ 



/ft \ 9/2 / i-t \ q/2 

El / \Yr\'^dr] +EI / \Zr\'^dr\ 



+ E 



\ g/2 



\f{r,0,0)\'dr\ 



+ \t- 



19/2 



sup EIZ^I" 

0<r<T 



<C\t-s 



q/2 



The proof is complete. □ 

Remark 2.8. From Theorem 2.6 we know that {{DgYt, DgZt)}o<e<t<T 
satisfies equation (2.21) and Zt = DtYt, fi x P a.e. Moreover, since (2.14) 
and (2.16) hold, we can apply the estimate (2.2) in Lemma 2.2 to the linear 
BSDE (2.21) and deduce supQ<«2'E|Zt|'? < 00. Therefore, by Lemma 2.7, 
the process Y satisfies the inequality (2.30). By Kolmogorov's continuity 
criterion this implies that Y has Holder continuous trajectories of order 7 
for any 7 < i - i. 

2.4. Examples. In this section we discuss three particular examples where 
Assumption 2.2 is satisfied. 

Example 2.9. Consider equation (1.1). Assume that: 

(a) f{t, z) : [0, T] X M X M ^ R is a deterministic function that has uni- 
formly bounded first- and second-order partial derivatives with respect to y 
and z, and f{t, 0, 0)^ dt < 00. 
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(b) The terminal value ^ is a multiple stochastic integral of the form 



(2.31) e= / g{tu...,tn)dWt,---dWt^, 



'[0,T] 

where n>2 is an integer and g{ti,...,tn) is a symmetric function in 
L2([0,r]"), such that 



sup / g{ti,...,tn~i,u)'^dti---dtn-i<oo, 

)<n<T J[0,T]"-i 

sup / g{ti,...,tn-2,U,vfdti---dtn-2<00, 

<u,v<T J \0,T]"-^ 



0<u,v<T J [0,TY 

and there exists a constant L > such that for any u,v & [0,T] 

/ \g{ti,...,tn-i,u) - g{ti, . . . ,tn-i,v)\'^ dti ■ ■ ■ dtn-i <L\u-v\. 

From (2.31), we know that 

Du£, = n i g{ti,...,tn-i,u)dWt^---dWt„_^. 

J[0,T]"-i 

The above assumption implies Assumption 2.2, and therefore, Z satisfies the 
Holder continuity property (2.24). 

Example 2.10. Let = Co([0, 1]) equipped with the Borel cr-field and 
Wiener measure. Then, is a Banach space with supremum norm || • ||oo, 
and Wt = uj{t) is the canonical Wiener process. Consider equation (1.1) on 
the interval [0, 1]. Assume that: 

(gl) /(t, y, z) : [0, 1] X M X R — > M is a deterministic function that has uni- 
formly bounded first- and second-order partial derivatives with respect to y 
and z, and /J /(t, 0, 0)^ dt < oo. 

(g2) ^ = (p{W), where — >■ M is twice Frechet differentiable, and the 
first- and second-order Frechet derivatives 6ip and satisfy 

\ip{u;)\ + \\6ip{u)\\ + ||5V(^)II < Ciexp{C2||a;|r^} 

for all u; € and some constants Ci > 0, C2 > and < r < 2, where || • || 
denotes the operator norm (total variation norm). 

(g3) If A denotes the signed measure on [0, 1] associated with 5(p, there 
exists a constant L > such that for all <9 <6' <1, 

K\x{{e,e'])\p <L\9-9'\p/^ 

for some p>2. 
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(2.32) 



It is easy to show that D^^ = X{{0, 1]) and DuDq^ = v{{0, 1] x (u, 1]), where v 
denotes the signed measure on [0, 1] x [0, 1] associated with From the 
above assumptions and Fernique's theorem, we can get Assumption 2.2, and 
therefore, the Holder continuity property (2.24) of Z. 

Example 2.11. Consider the following forward-backward stochastic 
differential equation (FBSDE for short): 

Xt = Xo+ [ b{r,Xr)dr+ [ a{r,Xr)dWr, 
Jo Jo 

Yt = V\^j^ X^drj+j^ f{r,Xr,Yr,Zr)dr- ZrdWr, 

where b,a, ip and / are deterministic functions, and Xq € M. 
We make the following assumptions: 

(hi) b and a has uniformly bounded first- and second-order partial deriva- 
tives with respect to x, and there is a constant L > 0, such that, for any 

s,te[0,T], xGM, 

\a{t,x) -a{s,x)\ <L|t-s|^/^. 
(h2) supo<i<r{|Ki,0)| + \ait,0)\} < oo. 

(h3) ip is twice differentiable, and there exist a constant C > and a pos- 
itive integer n such that 



X?dt 



+ 



T 



X^dt 



+ 



T 



X?dt 



<c{i+\\x\\^r, 



where ||x||oo = sup{|x(t)|, <t <T} for any x € C([0, T]). 

(h4) f(t, X, y, z) has continuous and uniformly bounded first- and second- 
order partial derivatives with respect to x, y and z and f{t, 0, 0, 0)^ dt < oo. 

Notice that in this example, ^{X) = ip{J^ X^ dt) is not necessarily globally 
Lipschitz in X , and the results of [16] cannot be applied directly. 

Under the above assumptions, (hi) and (h4), equation (2.32) has a unique 
solution triple {X,Y,Z), and we have the following classical results: for any 
real number r > 0, there exists a constant C > such that 

r/2 



E sup \Xt\^ < oo, 

0<t<T 



E\Xt-XX < C\t 



for any t,s £ [0,T]. For any fixed {y,z) € M x M, we have DQf{t,Xt,y,z) = 
dxf{t,Xt,y,z)DgXt. Then, under all the assumptions in this example, by 
Theorem 2.2.1 and Lemma 2.2.2 in [14] and the results listed above, we 
can verify Assumption 2.2. Therefore, Z has the Holder continuity prop- 
erty (2.24). 

Note that in the multidimensional case we do not require the matrix aa^ 
to be invertible. 
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3. An explicit scheme for BSDEs. In the remaining part of this paper, 
we let vr = {0 = to < ^1 < ■ • • < in = ^} be a partition of the interval [0, T] 
and |7r| = maxo<i<n-i l^i+i —ti\. Denote Aj = ti+i — ti,0 < i < n — 1. 

From (1.1), we know that, when t G 

(3.1) Yt = Yt^^, + f{r, Yr, Zr) dr - Zr dWr. 

Jt Jt 

Comparing with the numerical schemes for forward stochastic differential 
equations, we could introduce a numerical scheme of the form 



J ti 

t G [ti, ti^i),i = n — 1, n — 2, . . . , 0, 
where £ L^{n) is an approximation of the terminal condition ^. This leads 
to a backward recursive formula for the sequence {Y^.'^ , Z^^^}Q<i<n- In fact, 
once Y}'^ and Z}'^ are defined, then we can find Y}'^ by 

and {Zr'^}ti<r<ti+i is determined by the stochastic integral representation 



of the random variable 



Although {Z}'^}t^<r<u+i can be expressed explicitly by Clark-Ocone-Hauss- 
man formula, its computation is a hard problem in practice. On the other 
hand, there are difficulties in studying the convergence of the above scheme. 

An alternative scheme is introduced in [16], where the approximating pairs 
(y^jZ'^) are defined recursively by 

(3.2) y,- = y-^^ + / , y-^^ , E z; dr | f^^^ ) a, 

- J^^^' Z^dWr, te[ti,ti+i),i = n-l,n-2,...,0, 

where, by convention, E(^^ Z^ drlTt^^-^) = when i = n — 1. In [16] 
the following rate of convergence is proved for this approximation scheme, 
assuming that the terminal value ^ and the generator / are functionals of 
a forward diffusion associated with the BSDE, 

fT 

(3.3) max E|yt. -y^^l^ + E / \Zt - Z]^\'^ dt < K\Tr\. 

0<i<n ' ' Jq 
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The main result of this section is the fohowing, which on one hand im- 
proves the above rate of convergence, and on the other hand extends terminal 
value and generator / to more general situation. 

Theorem 3.1. Consider the approximation scheme (3.2). Let Assump- 
tion 2.2 he satisfied, and let the partition vr satisfy maxo<i<n-i A^/Aj+i < Li, 
where Li is a constant. Assume that a constant L2 > exists such that 

(3.4) |/(t2, y, z) - f{ti,y, z)\ < L2\t2 - 

for all ti,t2 G [0, T] and ?/, 2: G M. Then there are positive constants K and 5, 
independent of the partition it, such that, if \tt\ < 6, then 

(3.5) E sup iFt-yi^p + E / \Zt- Z]'\'^dt<K{\^^\+E\^-C\'^)■ 

o<t<T Jo 

Proof. In this proof, C > will denote a constant independent of the 
partition vr, which may vary from line to line. Inequality (2.24) in Theo- 
rem 2.6(b) yields the following estimate (Theorem 3.1 in [16]) with p = 2: 

J^E / " {\Zt-Ztf + \Zt-Zt,^,\^)dt<C\7r\. 

i=0 

Using this estimate and following the same argument as the proof of Theo- 
rem 5.3 in [16], we can obtain the following result: 

(3.6) maxE|yt,-y,7|2 + E ^ \Zt - Z^\'' dt < Ci\7T\ - C?)- 

0<J<n Jq 

Denote 

{0, if i = n; 
E Z:^dr\Tt,^, if i = n - 1, n - 2, . . . , 0. 

If ti <t < ti^i, i = n — l,n — 2, . . . ,0, then, by iteration, we have 

''ti+1 



(3.8) 



n ..T 

= C+ ^ fitk,Y,l,Zl)Ak-i- Z^dWr. 
1 — 11 t 



k=i+l 

Therefore, 



k=i+l 
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We rewrite the BSDE (1.1) as follows: 
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(3.9) 



where 



Yt=i+ [ f{r, Yr, Zr) dr- [ dWr 
Jt Jt 



k=i+l 



\Rt 



< 



/I 
fir,Yr,Zr)dr- J2 f{tk,Yt„Zt,)A,^ 
k=i+l 

" ftk ft 
^ / [f{r,Yr,Zr)-fitk,Yt„ZtJ]dr- f{r,Y,.,Zr)dr 

r \f{r,Yr,Zr)- fitk,Yt„Zt,)\dr+ T'"' \f{r,Yr,Zr)\dr. 

—,-1-1 •'*fe-l 



k=i+l 

By Lemma 2.2 and the Lipschitz condition on /, we have 

rT \ p/2 



eH \f{r,Yr,Zr)\^dr\ < oo, 



and hence, 



(3.10) 



E max 

0<i<n-l \ Jf 



k+i \ P 

\f{r,Yr,Zr)\dr 

i 

rT 



<\7t\p/'e(^I^ \f{r,Yr,Zr)\'dr^ 



\ p/2 



Define a function {t(r)}o<r<r by 

\ti+i, iiti<r <ti+i, z = n-l,...,0. 

By the Holder inequality, the boundedness of the first-order partial deriva- 
tives of /, (3.4), (2.24), Remark 2.8 and (3.10), it is easy to see that 



E sup \RJ\P<2P-^ 

0<t<T 



+ E max 

0<i<n-l V // 



ti+i \ P- 

\f{r,Yr,Zr)\dr 



(3.11) <(2r)f-iE r\f{r,Yr,Zr)-f{t{r),Yt^,),Zt^,))\Pdr 

Jo 
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/ rT \p/2 

<Cj7r|P/^ 

where, by convention, Rt = 0. In particular, we obtain 
(3.12) E sup \R1\'^ <C\-k\. 

0<t<T 

To simplify the notation we denote 

SY^'' = Yt-Yt'', 5Z^ = Zt-Z^ foralHG[0,r] 

and 

Zl = Zu-Zl fori = n,n-l,...,0. 
Then, when ti <t < ti+i, by (3.8) and (3.9) we can write 



fe-i 



k=i+l 



- 6Z^dWr + Rt +6C, 
where 6^'^ = ^ — ■ Therefore, we obtain 

(3.13) 6Y,^ = e( [fitk,Yt^,ZtJ-f{tk,Y,l,Zl)]Ak^, + RJ + 5C 

\k=i+l 



Denote = f{tk,Yt,, Zt,) - fitk,Y,l, Z^). From equality (3.13) for < t < 
tjj^i, where i < j < n — 1, and taking into account that 6Yf = 5Yf^ = 6(^'^ , 
we obtain 



sup \5Yt^\< sup E V \fX\^k-i+ sup \R;\ + \5C\ 

U<t<T U<t<T 0<r<T 



The above conditional expectation is a martingale if it is considered as a pro- 
cess indexed by t € [ti,T]. Thus, using Doob's maximal inequality, we obtain 

2 



E sup \6Y^''\'' <K sup 
U<t<T U<t<T 



E V \fl\Ak^i+ sup \R;\ + \6C\ 
<CE[ Y \ftj^k~i+ sup + 



\k=i+l 



0<r<T 
2 



<c\e[ Y +E^sup^|i?-|2+E|5e 



7r|2 



\k=i+l 
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From (3.12), we deduce 

E sup \SY,^\''<cIe( \ftj^k-i] + 



ti<t<T 



\k=i+l 



Using the Lipschitz condition on /, we obtain 

E sup \5Yi''\'^ <c\{T -ti)^E sup \6Ytl\'^ 

U<t<T [ i+l<k<n 

(n-l \ 2 

k=i+l J 

+c{E\sef + w\). 

Notice that 

n—l \ 2 / n—1 

A,, /. ^ 

Kk=i+1 / \k=i+l 

n-l 



/ "Zl A. , ftk+i 

E ^ / n\zt,-z:\\Tt,)du 

<L\e[ / E{\Zt,-Zl\\Tt,)du\ 

<2L?W |] l''^'E{\Zt,-Zum,)du 



(3.15) 



n— 1 



'2/ 



+^( E jy'n\Zu-zi\\Ft,)du 



= 2Li{h + h). 
Now the Minkowski and the Holder inequahties yield 

Izl ( ftk+i „ ^1/2 ^2 



/i<e( ^ I r^'(E(|Zi,-Zj|7-iJ)2dn| A 

<(r-t,) ^ r^'E(E(|Zi,-Zj|J-iJ)2ciu 

(3.16) ^ 

<iT-U) / E|Zi,-Z„ 



k=i+l 
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<C{T-ti) I \tk-u\du<C\'K\. 
In a similar way and by (3.6), we obtain 

n—l „f 



h<{T-ti) J2 / nzu-z:i\^du 
k=i+i 

(3-17) 

= {T-ti) 1 ^5Zl\^du<C\^\. 



On the other hand, 

(3.18) E(Z,';^A„„i)2 = E|Ztj2|A„_i|2<C|7r|2. 

From (3.14)-(3.18), we have 

E sup \5Yt''\^ <Ci{T -tif¥. sup l^y^^p 

ti<t<T i+l<k<n 

(3.19) 

+ C2(E|5ri' + |vr|), 

where Ci and C2 are two positive constants independent of the partition vr. 

We can find a constant 5 > independent of the partition vr, such that 
Ci(3(5)^ < I and T > 26. Denote I = [^] {[x] means the greatest integer no 
larger than x). Then / > 1 is an integer independent of the partition vr. If 
|vr| < 5, then for the partition vr we can choose n — 1 > ii > 12 > • ■ ■ > ii >0, 
such that, T-25e (tii_i,tij, T - 45 e {ti^-i^U^], . . . ,T - 251 € [0,tjj (with 
t-i = 0). 

For simplicity, we denote tj,, = T and = 0. Each interval [tij^-^ , ti.],j = 
0, 1, . . . , /, has length less than 36, that is, \ti. — ti._^^ \ < 36. On each interval 
[^ij+i ) iijjj j = 0, 1, . . . , /, we consider the recursive formula (3.2), and (3.19) 
becomes 

E sup \6Y,^\^<Ci{U^-U^^,fK sup \6Y,l\^ 

(3.20) 

+ C2(E|5F-|2^|vr|). 



Using (3.20), we can obtain inductively 



E sup \6Yt 



7712 



< Ciiu^ - t.,+j'E sup \6Yti\' + C2msyti r + ki) 
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+ C2(E|<5y,:,|2^|7r|) 



-K |2 
3 



< (Ci(3(^)2)*^-*^+iE|(5y7 p 
+ C2(E|5y-|2^|7r|) 

X (1 + Ci(3(5)2 + (Ci(3(^)=^)2 + • • • + (Ci(3<5)2)'^-'^+i) 

< E|,5y/ ? + 2C2(E|5y7 |2 + Ivrl) 
= (2C2 + l)E|<5yi" |2 + 2C72|7r|. 

By recurrence, we obtain 
E sup \SY^f 

< (2C2 + iy+^E\5C\^ + C2|7r|(l + (2C2 + 1) + • • • + (2C2 + 1)^') 

(3.22) 

< (2C2 + iy+^E\6C\^ + C2|7r|(l + (2C2 + 1) + • • • + (2C72 + 1)') 

Therefore, taking C = '^^^'^^^^^ — , we obtain 

E sup Idy^'^P < maxE sup l^^y^^p < C(|7r| + E|e - Tl^)- 

0<f<T 0<j<l ti.^,<t<ti. 

— — j+i— — 3 

Combining the above estimate with (3.6), we know that there exists a con- 
stant A' > independent of the partition vr, such that 

rT 



E sup lyt-yj'^p+E / \Zt- Z^\^dt<K{\7^\+E\^-^' 
o<t<T Jo 



7^|2^ 

□ 



Remark 3.2. The numerical scheme introduced before, as other similar 
schemes, involves the computation of conditional expectations with respect 
to the cj-field J~ti+i- To implement this scheme in practice we need to ap- 
proximate these conditional expectations. Some work has been done to solve 
this problem, and we refer the reader to the references [2, 4] and [8]. 

4. An implicit scheme for BSDEs. In this section, we propose an im- 
plicit numerical scheme for the BSDE (1.1). Define the approximating pairs 
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(Y'^,Z'^) recursively by 

_ en 

(4.1) = Y,l^ + /(^i,+i,y^-^^, i- Z^dr) A, - Z:^dWr, 

te [ti,ti+i),i = n- l,n-2,...,0, 

where the partition tt and Aj,i = n — 1,...,0, are defined in Section 3, and ^'^ 
is an approximation of the terminal value ^. In this recursive formula (4.1), 
on each subinterval [ij, tj+i), z = n— 1, . . . , 0, the nonlinear "generator" / con- 
tains the information of Z'^ on the same interval. In this sense, this formula is 
different from formula (3.2), and (4.1) is an equation for {{Y^^ , Z^)}ti<t<ti+x ■ 
When |7r| is sufficiently small, the existence and uniqueness of the solution 
to the above equation can be established. In fact, equation (4.1) is of the 
following form: 

(4.2) Yt = ^ + g(^j Zrdr^ - ZrdWr, te[a,b] andO<a<b<T. 

For the BSDE (4.2), we have the following theorem. 

Theorem 4.1. Let 0<a<b<T andp > 2. Let ^ be J^b-measurable and 
^ € L^(r2). Lf there exists a constant L > such that g:{^ x M, J7, (X" ^B) — )• 
(]R,,B) satisfies 

\g{zi) -5(^2)1 <L\zi - Z2\ 

for all zi,Z2 € K and g{0) € L'p{Q), then there is a constant 5{p,L) > 0, such 
that, when b — a< 5{p,L), equation (4-2) has a unique solution {Y,Z) € 
S^[a,b])xHma,b]). 

Proof. We shall use the fixed point theorem for the mapping from 
Hjr{[a,b]) into Hjr{[a,b]) which maps z to Z, where {Y,Z) is the solution of 
the following BSDE: 

(4.3) Yt=C + g(^j^ Zr dr^ - j^^ Z, dWr, te [a, b]. 

In fact, by the martingale representation theorem, there exist a progressively 
measurable process Z = {Zt}a<t<b such that E Zf dt < 00 and 

rb \ / / rb \ 

J-J + / ZtdWt. 



^ + 9\^J Zrdrj=E\^C + g\^J Zrdr 

By the integrability properties of S,,g{0) and z, one can show that Z G 
HP,{[a,b]). Define Yt = E{C + g{f^ Zrdr)\Tt),t £ [a,b]. Then {Y,Z) satisfies 
equation (4.3). Notice that y is a martingale. Then by the Lipschitz condi- 
tion on g, the integrability of £,,g{0) and z, and Doob's maximal inequality, 
we can prove that Y G S^{[a, b]). 
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Let be two elements in the Banach space H^{[a, 6]), and let {Y^^^ 
(y^,Z^) be the associated solutions, that is, 

Y;=^ + g( [ zl.dr\- f ZidWr, t e [a,b],i = 1,2. 



Denote 



Y = Y' -Y' 



Then 

(4.4) Yt = g(^l'zldr ]-g{ / zf dr 

for all t S [a, b] . So 



z = z^ — . 



ZrdWr 



Yt = K(g( / z^^dr] - g( / z^: dr 



Tt 



for all t S [a, 6]- Thus by Doob's maximal inequality, we have 



E sup \Ytf = sup 

a<t<h a<t<b 



(4.5) 



< CE 



< CE 



g(^j z}.dr^ ~^(/ 



dr 



dr — z^ dr 



<C{h-af'^^(^j^ \zrfdr^ 



where C > is a generic constant depending on L and p, which may vary 
from line to line. From (4.4), it is easy to see 

Yt = Ya+ [ ZrdWr 
J a 

for all t € [a, b] . Therefore, by the Burkholder-Davis-Gundy inequality 
and (4.5), we have 



fb _ \p/2 

E( / \Zr\^dr < CE sup 

/ a<t<b 



f Zr 

J a 



dWr 



(4.6) 



<C 



E|ya|P + E sup \Yt\P 

a<t<b 



nb \ p/2 

<C{b-ay'^E( / \zr\^dr] , 
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that is, 

\\Z\\Hv<Ci{h-af/'^\\z\\Hv, 

where Ci is a positive constant depending only on L and p. 

Take 5{p,L) = 1/Cf. It is obvious that the mapping is a contraction when 
b — a < 6{p, L), and hence there exists a unique solution (Y, Z) E Sjr{[a, b]) x 
Hma,b]) to the BSDE (4.2). □ 

Now we begin to study the convergence of the scheme (4.1). 

Theorem 4.2. Let Assumption 2.2 be satisfied, and let tt be any parti- 
tion. Assume that € U'[Vl) and there exists a constant Li > such that, 
for all ti,t2 G [0,T], 

\f{t2,y,z)- f{ti,y,z)\<Li\t2-ti\^'\ 

Then, there are two positive constants 5 and K independent of the parti- 
tion TT, such that, when |7r| < 5, we have 

/ rT \p/2 

sup |yi-y-|p + E / \Zt-z^\^dt] <K(|7r|p/2 + E|e-rr). 



E 

0<t<T 



Proof. If |7r| < 6{p,L), where S{p,L) is the constant in Theorem 4.1, 
then Theorem 4.1 guarantees the existence and uniqueness of (V^ , Z'^). De- 
note, for i = n — l,n — 2, . . . ,0, 

Zl , = / dr. 

Notice that {^j^, }i=n-i,n-2,....o here is different from that in Section 3. Then 

Z'^dWr, i = n-l,n-2,...,0. 
Recursively, we obtain 



fc=j+i 



Z^dWr, i = n-l,n-2,...,0. 



U 



Denote 

SC = ^-C, 6Y,^ = Yt-Y,\ 5Z^ = Zt-Z^, t€[0,T], 



and 



Zl = Zt^-Zl, i = n-l,...,0. 
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If t € i = n — 1, n — 2, . . . , 0, then by iteration, we have 
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(4.7) 



6Y^^ = 5C+ [fitk,Yt„ZtJ-f{tk,Y,l,Zl)]Ak- 

k=i+l 



5Z^.dWr + Rl, 



u 



where is exactly the same as that in Section 3. 

Denote = f{tk, Yt, ,ZtJ- f{tk, Ytl, Z^ ) . Then for t G [t^, f = n - 1, 
n — 2, . . . , 0, we have 

(4.8) 5Y,^ = e(6C+ Y1 fl^k-i + RlM. 

\ k=i+l I 

From equahty (4.8) for tj <t < tj+i, where i < j < n — 1, and taking into 
account that 6Y^ = 6Y^^ = 5^,'^, we obtain 

/ n \ 



sup \bYn< sup E V |/t'^.|Afe_i+ sup \K\^W\ 

ti<t<T U<t<T \f~^^ 0<r<T J 

The above conditional expectation is a martingale if it is considered as a pro- 
cess indexed by t for t G [ti,T]. Using Doob's maximal inequality, (3.11), and 
the Lipschitz condition on /, we have 

E sup \SY^\v 

ti<t<T 



<E sup 

ti<t<T 



E Yl i/rjAfc-i+ sup m + \6c\ 



<CE V |/rjAfc_i+ sup \K\ + \SC\ 



<c|e|^^ +E^sup^|i?-r + E|5rr| 



\k=i+l 



\k=i+l 



<C{{T-UrE sup \6Y,l\ 



i+l<k<n 



+E Y \^tj^k-i] +\7r\p/^+n6c\n, 

\k=i+l J ) 
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where, and in the fohowing, C > denotes a generic constant independent 
of the partition vr and may vary from line to Une. On the other hand, we 
have, by the Holder continuity of Z given by (2.24), 



\k=i+l 



\k=i+l 



z 



1 







r Z^dr 




ltk-1 





" rt 



\k=i+l ''^k-l 



<A E / \Zt,-Zr\dr+ 



I Z'p — Zy. I dv 



< C7|7r|P/2 + 2P"^e(^£ \Zr - Z^ldr^ 

/ rT \p/2 

< C\tt\p/^ + 2P-\T - ti)P/^Ei^J^ \Zr - drj 

/ rT \p/2 

= C\TT\P/^ + 2P~\T-U)P/^Ei^J^ \5Z;\^drj . 



Hence, we obtain 



(4.9) 



E sup ijy^^ip 

U<t<T 

<cJiT-tirE sup \6YtJ 

L i+l<fc<n 



/ rT \p/2 

+ (T-t,)P/2E^ \5Z^\^drj 

+ |7r|P/2+E|5rr|, 



where Ci is a constant independent of the partition vr. By the Burkholder- 
Davis-Gundy inequality, we have 

p/2 



(4.10) E 
From (4.7), we obtain 



\5Z:Vdr\ <CpE 



6Z^ dWr 



(4.11) r 5Z^^dWr = 5i^+ flAk^i + Rl-5Y,^ 

k=i+l 
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Thus, from (4.10) and (4.11), we obtain 

rT \ p/2 

') 



E 



rT \ ) 

\6Z^\^dr\ 



k=i+l 



+E\6c\^ +^Rt\'' +^sYt''\p y 



Similar to (4.9), we have 

/ f-T \p/2 



<C2{{T-tiYE sup \5YtJ 

i+l<k<n 



where C2 is a constant independent of the partition vr. 
If C2(T - ti)P/^ < i, then we have 



E 



rT \ p/2 

/ \6Z;\^dr) <2C2{T-tiYE sup \5YtJ 

Jti J i+l<k<n 



(4.12) 

+ 2C2i\7T\P/^ + E\5Cn- 
Substituting (4.12) into (4.9), we have 

E sup \6Yt''\P 

ti<t<T 

<Ci{l + 2C2{T-UY'^){T-tiYE sup \5Yt 



(4.13) 



i+l<k<n 

+ Ci{l + 2C2{T - tiY^^){\TT\P/^ + E\6C\n 
<2Ci{T-UYE sup \5YtJ + 2C^{\T:\P/^+E\5en. 

i+l<k<n 



We can find a positive constant 5 < 6{p, L) independent of the partition vr, 
such that. 



C2(35)P/2 



< 



2' 



2Ci(3(5)P < i 



(4.14) 

(4.15) ^2 

and T > 25. Denote I = [^]. Then / > 1 is an integer independent of the 
partition vr. If [vrl < 5, then for the partition vr we can choose n — l>ii>i2> 
■■■>ii>0, such that, T-25e (ti,_i,tij, T - A6 e {ti^-i,ti^], . . . ,T - 251 G 
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[0,tjj (with t^i = 0). For simplicity, we denote = T and ij^^j = 0. Each 
interval [ti■_^^-^,ti.],j = 0,1, . . . ,l,haslength. less than 3(5, that is, {U- —ti_^^-^\ < 
36. On [ti._^j^,ti.], we consider the recursive formula (4.1). Then (4.13)-(4.15) 
yield 

E sup \6Yt''\P 

< 2Ci{U^ - U^^.YE sup \6YtJ + 2Ci(Kr/2 +E|5y- |f) 

ij^l+l<k<ij 

(4.16) 

< 2Ci{3Sf-E sup \5Yt^ \p + 2Ci(|7r|P/2 + K\SYtl \p) 

<\ sup \5YtJ + 2C,{\i:f^+n5Y,l\n. 
^ jj+i+i<fc<jj 

As in the proof of (3.21) and (3.22), we have 

E sup |(5yt^|*'<(4C7i + l)E|5Ft7 |P + 4Ci|7r|P/2 

and 



j+l ] 



E sup |5y-|p<^il^il±ll^(E|5ri' + M^/') 

Therefore, we obtain 

E sup \SY^f <mscKY. sup l^yj'^IP 

(4.17) 



Q<t<T 0<j<l U^^^<t<U^ 



E 



On [tj^.^-^, tj^.], J = 0, 1, . . . , Z, based on the recursive formula (4.1) and (4.17), 
inequality (4.12) becomes 

' rU^ \p/2 
\5Z;\^drj 

<2C2(ti,. -ti^.+J^E sup \5YtJ + 2C2i\7r\P/^+n6Cn 

^j-^-l+l<k<ij 

<2C2{3S)PE sup \SYtJ + 2C2{\7T\P/^+E\5C\n 

ij_l_l + l<fc<ij 

<^E sup \6YtJ + 2C2{\7T\P/^+n6Cn 

< (^^^^i±^+2C72)(i7rr/2+E|5rn. 
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Thus 



E 



\ P/2 

\5zn'dtj 



(4.18) 



\j=0 J 

<(/ + 1)p/2-i^]e(' I''^ \5Z^\^dt\ 



'J 

< (/ + i)p/2 |^ 3(4Ci + 1)^+1 ^ 2^^^! ^ 
Combining (4.17) and (4.18), we know that there exists a constant 

independent of the partition vr, such that 

cT \ p/2 

, ^ , O - \ 



sup |yt-yt"|p + E / \Zt-zj\^dt\ 

<t<T \Jo J 

<K{\ir\P/'^ +E\C-Cn- □ 



E 

0<t<T 



Remark 4.3. The advantages of this imphcit numerical scheme are: 

(i) we can obtain the rate of convergence in sense; 

(ii) the partition vr can be arbitrary (|7r| should be small enough) without 
assuming maxo<i<n-i Ai/A,j+i < Li. 

5. A new discrete scheme. For all the numerical schemes considered in 
Sections 3 and 4, one needs to evaluate processes {Z^}Q<t<T with continuous 
index t. In this section, we use the representation of Z in terms of the 
Malliavin derivative of Y to derive a completely discrete scheme. 

From (2.21), {DgYt}o<0<t<T can be represented as 

(5.1) DeYt = E(^pt,TDei + pt,,L»e/(r,y,, Z,) dr| J^j , 
where 

(5.2) pt^r = expj^ Ps dWs + ^ (^as- ^^^^ ds^ 
with as = dyf{s,Ys,Zs) and /3s = dzf{s,Ys, Zg). 
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Using that Zt = DiXt^ jJ- x P a.e., from (1.1), (5.1) and (5.2), we propose 
the following numerical scheme. We define recursively 



YZ = nYZ,, + fik^i,Y,l, , ZlJA,\TtJ, 



\ k=i / 

i = n — l,n — 2, . . . ,0, 

where pj^ = l,i = 0,1, . . . , n, and for < i < j < n. 



(5.4) 



. k=i 



An alternative expression for Pt^^tj given by the following formula: 
Pl,t,=^^AY.^.f{tkXtl,Zl){Wt,^^-Wt,) 



(5.5) 



+ E ( /(^'^^ - l[dJitk,YZ,Zl)]A A 



However, we will only consider the scheme (5.3) with p'^_ ^. given by (5.4). 
We make the following assumptions: 

(Gl) f{t,y,z) is deterministic, which implies DQf{t,y,z) = 0. 
(G2) f{t,y,z) is linear with respect to y and z; namely, there are three 
functions g{t), h{t) and /i(t) such that 

f{t,y,z)=g{t)y + h{t)z + fi{t). 

Assume that g, h are bounded and /i G L^([0,T]). Moreover, there exists 
a constant L2 > 0, such that, for all ti,t2 £ [0,T], 

\git2) - g{ti)\ + \h{t2) - hiti)\ + 1/1 (t2) - fiih)\ < L\t2-h\'/\ 

(03) Esupo<e<T l^eCr < for ah r > 1. 
Notice that (Gl) and (G2) imply (ii) and (iii) in Assumption 2.2. 

Remark 5.1. We propose condition (Gl) in order to simplify 
{Z^.}i=n-i,...,o in formula (5.3). In fact, there are some difficulties in general- 
izing the condition (G)s, especially (Gl), to a forward-backward stochastic 
differential equation (FBSDE, for short) case. 
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If we consider a FBSDE 

Xt = XQ+ I b{r, Xr) dr + I a(r, X,) dWr, 
^ Jo Jo 

Yt = ^+ I f{r,Xr,Yr,Zr)dr- I ZrdWr, 
^ Jt Jt 

where Xq G M, and the functions 6, u, / are deterministic, then under some 
appropriate conditions [e.g., (hl)-(h4) in Example 2.11] Z^. for i = n — 
1, . . . , in (5.3) is of the form 

Zl=^\pl_^^,^Dui 



n-l 



+ p7 f d^f(tk+i,x7 ,y,'^ ,zf )Dt-xr a* 



k=i 



U 1 



where {X'^ ^V" , Z^) is a certain numerical scheme for (X, Y, Z). It is hard to 
guarantee the existence and the convergence of Malliavin derivative of X^ , 
and therefore, the convergence of Z'^ is difficult to derive. 

Theorem 5.2. Let Assumption 2.2{i) and assumptions (G1)-(G3) be 
satisfied. Then there are positive constants K and 6 independent of the par- 
tition IT, such that, when j-Trj < 5 we have 

/ 1 \P/2 

E max {\Yt^ - Y^J + \Zt^ - Z^f} < i^|^|P/2-p/(2iog(i/W)) / j^g _\ ^ 

Proof. In the proof, C > will denote a constant independent of the 
partition vr, which may vary from line to line. Under the assumption (Gl), 
we can see that 

(5.6) Zl = E{pl^^^,^Dt^^\TtJ, i = n - 1, n - 2, . . . , 0. 

Denote, for z = n — 1, n — 2, . . . , 0, 



dZ]'. = Zt. - Z 



til 



6Y,^^=Yt^-Y,^. 



Since |e^ — e^| < (e^ + ey)\x — y\, we deduce, for all i = n — 1, n — 2, . . . , 0, 

<n\pu,tn-pi+„tj\DtA\\j'u) 



<E(^\DtMpu,tn+pi^,,t:, 



1 



h{r)dWr+ / g{r)dr-- I h{rY dr 



T 
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/ h{r)dWr- / 9{r)dr 

k=i+l •'tfc 



k=i+l 



1 f«fc+i \ 



/i(r) dWr 



+ 



|5r(r)|dr 



+ 2 



h{rfdr 



From (G2), we have 



n-l 



fc=i+l 



tfe+i 



hirfdr 



<CA sup iL'sCOf «^P| / /^(O'^W^r jY 

^o<6i<T ^ Vo<t<T yh )J 

where Ci > is a constant independent of the partition vr. 
In the same way, we obtain 

\Dui\Pu,t^<Ci{ sup \Dei\)( SUV expj C h{r)dwX]. 

^o<e<T \o<t<T Ut }/ 



Thus for i = n — 1, n — 2, . . . , 0, 



\5ZI\<2Cie(( sup ll^eCl)!^ sup expj / h{r)dW, 
\^o<e<T ' \o<t<T [Jt 



h{r) dWr 



+ 



ti+i 



\g{r) \ dr -\ — / h{r) dr 



<2Cx^{{ sup \D6i^{ sup expj / h{r)dwX\ 
W<e<T \o<t<T [Jt jj 



sup 

L0<fc<n-1 



h{r) dWr 



+ sup 

0<A;<n-l Jtk 



tk+l 



\g{r)\dr 



+ n sup 



0<fc<r!,-l Jtj. 



•fe+i 



/i(r)^ dr 



MALLIAVIN CALCULUS, NUMERICAL SOLUTION OF BSDE 



41 



The right-hand side of the above inequahty is a martingale as a process 
indexed by i = n — 1, n — 2, . . . , 0. 

Let r]t = exp{— h{u) dWu}- Then, r]t satisfies the fohowing hnear stochas- 
tic differential equation: 

di]t = -h{t)rit dWt + \h{tfr]tdt, 
% = 1- 

By (Gl), (G2), the Holder inequality and Lemma 2.4, it is easy to show 
that, for any r > 0, 

cT 



e( 



sup exp 

Vo<i<T 



h{u) dWu 



:E(^exp| [ h{u)dWu\ sup exp|- [ h{u)dWu 



(5.7) 



< Eexp 



h{u) dWu 



{-I 

X (e sup exp<^ -2r / h{u)dWu } 
V o<t<T I Jo J 



1/2 



h{u) drUE sup rjf) < oo. 
J ^ o<t<T ' 

For any p' € (Pi f)' Doob's maximal inequality and the Holder inequality, 
(G3) and (5.7), we have 

E sup \bZl\^ 

0<i<n 



<CE( sup |Z?e^|) ( sup exp<{ / h{r)dWr 



0<t<T 

k + 1 



sup 

0<fc<r!.-l 



sup 

0<fc<n-l Jt 



h{r) dWr 



tk+1 I 

\g{r)\dr + - sup 



0<fc<n-l Jtfe 



IP 



h{r)'^ dr 



<C 



E( ( sup \DeC\ 

0<9<T 



PP'/(p'-p) 



X ( sup exp<j / h{r)dWr 

vO<t<T 



pp'/{p'~p) 







/■<fc+i 




X 


El sup 


/ h{r) dWr 


+ sup / 




\0<fc<n-l 


Jtk 


0<k<n-lJtk 



ip'-p)/p' 



\g{r)\ dr 
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<C 



E( sup \Dei\ 
o<e<T 



2pp'/{p'-p) 



T 



E(^supexp|y^ h{r)dWr^j 



1 rtk+i \ P 
H — sup / /i(r)^ dr 

2 0<fc<n-l Jtfe 
p'/{2{p'-p)) 

N 2pp7(p'-p)-|p'/(2(p'-p)) 



p/p' 







ftk+i 


i 


X 


E sup 


/ /i(r) dWr 


+ E sup 




0<fc<n-l 


Jtk 


0<fc<71-l V 



tk+l \ p 

\g{r)\dr 



tk+l 



+ E sup 

0</fc<n-l \Jtk 

For any r > 1, by the Holder inequality we can obtain 



/i(r)^ dr 



p/p' 





ftk+i 


P' r 


ftk+i 




/i = E sup 


/ /i(r) dWr 


<<E sup 


/ h{r) dWr 


p 


0<A,<n-l 


Jtk 


I 0<fe<n.-l 


Jtk 





< 



. fc=0 



tk+l 



h{r) dWr 



tk 



p r 



1/r 



For any centered Gaussian variable X, and any 7 > 1, we know that 

E|Xp<C'V/^(E|^P^/^ 
where C is a constant independent of 7. Thus, we can see that 

l/r 



Take r ■ 



21og(l/k|) 



Assume |7r| is small enough; then we have 



p72 



It is easy to see that 



sup 



0<fc<n-l \Jtk 



tk+l \p 

\g{r)\dr] <C\7r\P 



and 



/3 = E sup 

0<fc<n-l \Jtk 



■k+l 



h{r)'^dr] <C\tt\- 
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Consequently, we obtain 

(5.8) E sup <C7|7r|P/2-p/(2iog{i/kl)) 1 y 

0<j<n ' V FI/ 

Applying recursively the scheme given by (5.3), we obtain 
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p/2 



\ k=i+l 

Therefore, for i = n — l,n — 2, . . . ,0, 



i = n — 1, n — 2, . . . , 0. 



\fc=i+l 

where i?^ is exactly the same as in Section 3 and = C ~ C = 0- In fact, 
we keep the term bi^ to indicate the role it plays as the terminal value. 
For J = ?i — 1, n — 2, . . . , i, we have 



|5y,^|<E E l/(ifc,^t.,^tJ-/(tfe,^i:,^t:)|Afc-i 

\fe=i+l 

+ sup \Rl\-rW\ tA. 

0<t<T J 
By Doob's maximal inequality and (5.8), we obtain 
E sup \bYt''.\P 

i<j<n 

<Ce( E \fitk,Yt„Zt^) - f{t,,Y,l,Zl)\A,^ 

\k=i+l / 

<C E E \yt,-Ytl\^u-i] +E E \Zt,-Zl\Ak^-, 



Kk=i+l 



\k=i+l 



+c(|7r|p/2+E|5e"r) 

<C2(r-ti)^E sup l^t, 

j+l<fe<n 

+ C3(|7rr/2-p/(2iog(i/kl))(^log^ 



p/2 



+ E|(^^^| 



where C2 and C3 are constants independent of the partition vr. 

We can obtain the estimate for Emaxo<j<n|yj. — by using simi- 

lar arguments to analyze (4.13) in Theorem 4.2 to get the estimate for 

Esupo<t<T\yt-yt''\- □ 
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